Imaginary deuteron optical potential due to elastic and inelastic breakup 
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The contributions to the reaction cross section from the elastic and inelastic breakup processes, 
calculated within the post-form distorted wave Born-approximation theory, are used as constraints 
to determine the contributions to the imaginary part of the deuteron optical potentials (IPDOP) 
due to the breakup channels. The Coulomb part of this potential due to the elastic breakup process 
is seen to account for the long range absorption in the optical potential. The nuclear parts of 
the IPDOP due to the elastic and inelastic breakup modes peak in different regions of the nuclear 
surface, with the latter being almost an order of magnitude larger than the former. This makes 
the IPDOP due to the breakup channels determined by us stronger than those calculated earlier 
ignoring the inelastic breakup mode. 
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I. INTRODUCTION 

In collisions between two nuclei the breakup of the pro- 
jectile into two or more fragments is often a strong reac- 
tion channel, which affects not only the imaginary part 
but also the real part of the corresponding optical po- 
tential. This leads to a dynamical polarization potential 
(DPP) which has to be added to the real potential cal- 
culated by the double folding (DFM) models (see eg. a 
recent review Otherwise, the real part of the DFM 

potentials for weakly bound projectiles (eg. ^'''Li and 
^Be) require arbitrary re-normalization factors in order 
to fit their elastic scattering data ||^. Whereas, fold- 
ing model calculations have been performed for both real 
as well as imaginary parts of the nucleon optical poten- 
tials [|j , for the case of light ions they are confined only to 
the real potentials which together with a phcnomcnolog- 
ical imaginary part is used to describe the corresponding 
elastic scattering. 

One of the problems associated with the microscopic 
calculations of the imaginary part of the light ion optical 
potential has been to include the effects due to breakup 
of the projectile in the field of the target nucleus, which is 
a strong reaction channel for these nuclei. Experimental 
studies have shown that even for strongly bound projec- 
tiles the probability of breakup increases drastically with 
increasing beam energy For example, the cross 

section for breakup of the a particle into n -I- '^He in- 
creases by, at least, an order of magnitude as the beam 
energy is varied from 65 MeV to 140 MeV §,§. Thus 
the effects of breakup are important also for the tightly 
bound projectiles for beam energies above 30 MeV/A. 

The optical potentials due to the breakup channels 
have been calculated by several authors in the past 



-|11|. Most of them are based on the coupled chan- 
nels (CC) techniques where the excitation of the breakup 
channel and its feedback on the elastic channel is stud- 
ied. However, such calculations are rather complicated 
as one has to find reliable approximations to include the 
higher order effects and the complete breakup contin- 
uum in the calculations |l0|JTl| ]. Moreover, the inelastic 
breakup mode, which dominates the total breakup cross 
sections [|l^ can not be included in these calculations. 

In this paper we follow a method introduced in [p^ , 
where it was shown that unitarity of the scattering ma- 
trix makes it possible to investigate the influence of the 
breakup process on the elastic scattering even without 
introducing the coupling of the breakup channel back to 
the elastic channel. In this procedure, the elastic scatter- 
ing and breakup reaction are investigated separately. In 
the first step, the breakup of the projectile in the nuclear 
and Coulomb fields of the target nucleus is calculated fol- 
lowing the post form distorted- wave Born-approximation 
(DWBA) theory. In this first order theory, which repro- 
duces the experimental breakup data rather well, only 
the coupling of the elastic channel to the breakup channel 
is considered. The contribution of each partial wave of 
the incident projectile to the total breakup cross section 
can be explicitly determined within this theory. Without 
such a partial wave decomposition, the present approach 
would have not been feasible. 

In the second step, the elastic scattering of the projec- 
tile is calculated from the known optical potential. We 
determine the reaction cross section for each partial wave 
(which are uniquely determined by the imaginary part of 
the corresponding phase-shifts) and split it (using the 
unitarity of the scattering matrix) into two parts, one 
due to the breakup channels and the another due to the 
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rest. Since, the reaction cross sections out of a specific 
channel can be related to the expectation value of the 
imaginary part of the optical potential associated with 
that channel (which is calculated with the corresponding 
optical model wave function in the entrance channel) p4[ , 
we use the breakup cross sections calculated within the 
post form DWBA theory as constraints in a fitting pro- 
cedure to determine the imaginary part of the potential 
due to the breakup channels. We prefer not to call it 
as the dynamical polarization potential as this phrase is 
used for potentials having both real as well as imaginary 
parts. 

The formalism used in our calculations are discussed 
in the next section. The results and their discussions are 
presented in section III. The summary and conclusions 
of our work are given in section IV. 



II. METHOD OF CALCULATIONS 

We write the phenomenologically determined optical 
potential U{r) (= V{r) + iW{r)) as 



Uir) = [Uir)-Ubuir)]+UUr), 



(2.1) 



where Ubu is the dynamical polarization potential due 
to the breakup channels. The wave functions yi{r) and 
Wi{r), corresponding to potentials U{r), and [U{r) — 
Ubuii")] (= Ubareif), the bare potential) respectively, sat- 
isfy the following radial Schrodinger equations 

'Pyi{r) ^(^+l) i , . . 
— + [k -U{r) ] yi{r) = 0, (2.2) 



(Pvijr) 



+ [k^ - Ubareir) 



£(£ + 1), 



Wi{r) = 0, (2.3) 



where k is the wave-number of the incident deuteron. 
For r > Rp (where Rp is the distance beyond which the 
nuclear interactions can be ignored), the wave functions 
yi{r) and wi{r) are normalized according to 

yeir) ^ e'^' [cos SeF^ikr) + sin 6eGe{kr)], (2.4) 
we{r) - e'^° [cos S°fFe{kr) + sin S"iGi{kr)], (2.5) 

where Fi and Gi are the regular and irregular Coulomb 
functions. Si and are the scattering phase-shifts cor- 
responding to potentials U{r) and Ubare{r)- 

The expressions for the partial wave amplitudes can be 
written in either of the following two forms fl^ ] , 

1 f°° 

fi = -T Fe{kr)Ubare{r)yi{r)dr 
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ff + fF, 



Ft{kr)Ubu{r)vi{r)dr 



(2.6) 



and 



h = -T 



Fi{kr)Ubare{r)we{r)dr 



k Jo 



wi{kr)Ubuir)yi{r)dr 



(2.7) 



These decompositions of the partial wave amplitudes 
were used in an earlier study of the optical poten- 
tial due to breakup channels. However, the breakup am- 
plitude was assumed to be obtained from the difference 
fi — instead of fg — This neglected the fact that 
the amplitude from the bare potential is affected by the 
presence of the breakup potential and generated a de- 
pendence on the real part of the breakup potential. As 
will be shown below, a separation of the reaction cross 
section into contribution of various channels requires the 
knowledge of the real potential V{r) only. 

Although the reaction cross section (which includes 
contributions from all the inelastic channels) can be cal- 
culated directly from the partial wave amplitudes de- 
scribed above, we use here an expression where it is writ- 
ten in terms of the imaginary part of the optic al pot en- 
tial (110. This method was used earher ||l|,|l^g to 
calculate the contributions to the imaginary part of the 
optical potential from the fusion channels. We can write 
the reaction cross section (cta) as. 



2n 

hv 



<X^^\W\x. 



(+) 



(2., 



where v is the relative velocity in the entrance channel 
and Xi^^ the full solution of the Schrodinger equation 
(whose radial part is yi). We can also write 



l)T,. 



(2.9) 



where the transmission coefficient (Ti) is given by 



4 f°° 
T- / \ye{r 



)[^W{r)dr 



(2.10) 



It may be noted that can also be related to the am- 
plitudes fi (Eqs. (6)^ and (7)) by = 4(|/,|2 - //), 
where // denotes the imaginary part of fi. However, the 
advantage of Eq. (10) lies in the fact that it involves a 
linear dependence of Ti on the imaginary potential. This 
allows us to split into terms corresponding to the con- 
tributions from different channels, as will be discussed 
below. 

Using unitarity of the S'-matrix, the transmission co- 
efficient Tf can be written as 



^l~[Su[- 



(2.11) 



where S represents the scattering matrix and £ denotes 
the elastic channel. For simplicity of notation we take 
the projectile and target nuclei to be spin less; hence £ 



2 



corresponds to the total spin, and c describes any other 
channel with total angular momentum i. Thus Eq. (11) 
enables us to express the transmission coefficient Tg as 
a sum (or integral for continuous channels) over all the 
reaction channels. This allows us to write 



^bare 



(2.12) 



for each partial wave t. In this equation <Jh-up,d repre- 
sents the contribution to the reaction cross section from 
the breakup channels, while cr^"'''^ is the reaction cross 
section corresponding to the remaining channels. Fol- 
lowing Refs. p^Jl^, we decompose the total imaginary 
potential M^(r)into a bare component and a component 
due to breakup as, W(r) = Whareir) + VFfc„(r). Then, 
expressions similar to Eq. (9) can be written for a^^"^^ 
and (Tb-up.d with corresponding T^'s given by. 



and 



rpbu 

^ hv 



\yi{r)\^Wtu{r)dr, 



(2.13) 



(2.14) 



where the Wbu consists of a part due to the elastic 
breakup Wdiss and a part due to the inelastic breakup 
Winbu- In our fitting procedure, potentials Wdiss and 
Winbu (with a certain a priori assumed form) are var- 
ied so that the elastic or inelastic breakup cross sections 
(calculated within the post form DWBA theory which is 
described below) are reproduced for each partial wave. 
We impose the constraint that Wbu < W for all r. Of 
course, the potentials due to breakup so determined, are 
specific to our breakup cross section and it could be dif- 
ferent from such potentials defined by other authors. 

It may be noted that dependence on the real poten- 
tial entering into the transmission coefficients (Eqs. (10), 
(13) and (14)) is only through that of the phenomenolog- 
ical optical potential that is used to calculate yi. In the 
calculations presented in this paper we assume that this 
potential is known from the description of the elastic scat- 
tering. Thus no information about the real part of the 
potential due to the breakup channels can be extracted. 
However, had one started from a real potential calcu- 
lated within a double folding model, it would have been 
necessary to include a dynamical polarization potential 
(having a real part) in order to reproduce both the elastic 
scattering and breakup probabilities simultaneously. It is 
also worthwhile to note that if the energy dependence of 
the imaginary potential is known over a sufficiently large 
range of energies, the dispersion relations may be helpful 
in getting the corresponding real potential [^,^ . 

In the post form DWBA theory of the inclusive 
breakup reaction (eg. d + A ^ p + X, to be represented 
as {d,p)) the total breakup cross section is defined by 

0m 



^b-—up{d.p) 



(Pa{d, p) 



(2.15) 



where 



d^<j{d,p) 



is the double differential cross section for 



dVLpdEp 

the reaction {d^p), which is the sum of the elastic and 
inelastic breakup modes. The former (where X corre- 
sponds to n + A(g.s.)) is given by 



d^(j{elastic) 
d^pdEp 



p{phase) ^ I 13, 



(2.16) 



Using a zero range approximation, the amplitude Pe„m^ 
can be written as 

Pi^ra^^Do J dVx(,-)*(kp,-^r)F,„(fc„,r) 

X r;„™„(r)x^+Hkd,r)A(r)P(r). (2.17) 

Dq is the zero range constant for the d p+n vertex. Its 
value has been taken to be 125 MeV fm^/^ which is con- 
sistent with the known properties of the p-n system. The 
function A(r) takes into account the finite range effects 
within the local energy approximation (LEA) |p2|,^ . We 
have used the form of A(r) as that given in with a 
finite range correction parameter of 0.621, a value used 
in most of the calculations on deuteron induced transfer 
and breakup reactions. P{r) accounts for the nonlocality 
of the optical potentials. This is calculated by following 
the method of Perey and Buck |^,^ , with the nonlocal- 
ity parameters of 0.85 in the neutron and proton channels 
and 0.54 in the deuteron channel. 

In Eq. (17), x"*^ are the optical model wave functions in 
the respective channels with fc's being the corresponding 
wave numbers. Fi^{kn,r) is the radial part of the wave 
function in the n + A(g.s.) channel. In Eq. (16), p{phase) 
is the three-body phase-space factor p5| , p6[ . It should 
be noted that the integrand in Eq. (17) involves three 
scattering wave functions which are asymptotically oscil- 
latory. This makes the radial integrals involved therein 
very slowly converging. However, integrals of this kind 
can be effectivelyevaluated by using a contour integra- 
tion method [ P7y2q| . 

To calculated the cross section for the inelastic breakup 
process, where X can be any two-body channel of the 
B = n + A system, we start from a T matrix, 

'^Bxxi'^lVnplMdxld'^ >, (2.18) 

where (I)a and cfid denote the ground state wave functions 
of the target nucleus A and the projectile (deuteron) re- 
spectively, represents the complete scattering state 
of the system B with the boundary condition X. The in- 
tegration over the internal coordinates of in Eq. (18) 
leads to a form factor for the inelastic process. The calcu- 
lation of the form factor simplifies greatly if we use a sur- 
face approximation |p9[| , where we assume that the main 
contribution to Td.px comes from the region outside the 



Td,px = < '^ehp 
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range of the nuclear interaction. The validity of this ap- 
proximation has been tested by Kasano and Ichimura |3C|] 
by evaluating this integral without recourse to this ap- 
proximation. These authors find that the surface ap- 
proximation is valid for the deutcron induced breakup 
reaction even at the lower beam energy of 25 MeV. Thus 
we can represent the radial part of the form factor (F^) 
by its asymptotic form 



X {St 



where ji^ and h^^^ denote the spherical Bessel and Han- 
kel functions respectively. Si^i^ are the scattering matrix 
elements for the elastic channel corresponding to the an- 
gular momentum £„. Now, it is straight forward to carry 
out the integrations over the angles of the unobserved 
particle to get the double differential cross section from 
the triple differential cross sections. This leads to a re- 
duced T matrix for the process d + A ^ p + X , 



X 



2 V mxkx 
5t^xW\Kr), 



(2.19) 



d,pX 



rrinkn Sex 



mxkx Se^e„ - 1 

xDoJ dVx(-)*(kp,^r) 

X [F,„(fc„,r)-j,„(/c„r)]y£„Jf) 

X xi+'(kd,r)A(r)P(r). (2.20) 

In order to calculate the double differential cross section 
for the inelastic breakup, one has to sum over all the 
channels X ^ In- Since, in Eq.(20), the entire depen- 
dence on channel X rests solely in the S matrix S^^Xi 
this summation can be easily carried out using the uni- 
tarity of the S matrix 



.x\ 



^ 1 



15, 



(2.21) 



Therefore, we only need to know the S matrix elements 
of the elastic scattering to determine the double differ- 
ential cross sections for the inelastic breakup, which can 
be written as 



d^ (J (inelastic) 
d^lpdEp 



^f^h^ r.^\ \ ^ ( ^reaction i elastics 

plphase) (cr^^ /a^^ ) 
X I /3£„m„ - Plm^ f (2.22) 



In this equation o-^e^ction ^^^^ ^elastic -j-j^g reaction and 
elastic scattering cross sections for the neutron-target 
system corresponding to the partial wave £„ respectively. 
Pg m defined in the same way as Eq. (17) with the 
wave function Fg^ being replaced by the spherical Bessel 
function. More details of the derivation of the inelastic 
breakup cross sections can be found in Refs. | |T^ , |2^ 

The total cross section for the reaction (d,p) can also 
be written as 



^b—up(d,p) 



h—up{d,p) 



(2.23) 



where ^Pi'^'P) jg |-j^g transmission coefficient for the 
(d,p) breakup reaction, which is also termed as the 
breakup probability in ||l^. The total breakup prob- 
ability is given by 



rpb—up,d rpb~up{d,pn) 



rpb up(d,p) (■jj^g^Qg^jg^ 



_^rpb "P('*'")(^j^g;Qg^jg'j 



(2.24) 



In Eq. (24), T; 



b—up{d,pn) 



, represents the breakup prob- 



ability for the elastic breakup mode as defined above. 
The total breakup cross section <Jb-up,d is obtained from 



b—up,d 



by following an expression similar to Eq. (9). 



III. RESULTS AND DISCUSSION 

Apart from the zero range constant, finite range and 
nonlocality parameters described already in section II, 
we require the optical potentials in the deuteron, proton 
and neutron channels to calculate the breakup cross sec- 
tions. These have been taken from the global sets given 
by Daehnick, Childs and Vrcelj |^ (for the deuteron 
channel) and Becchetti and Greenlees ||3^] (for the pro- 
ton and neutron channels) respectively. 

In Fig. 1 we show the results for the breakup prob- 
ability for the deuteron incident on a ^^V target at the 
beam energy of 56 MeV, calculated within the post form 
DWBA theory. In this figure we have also shown the 
total transmission coefficients calculated with the same 
deuteron optical potential. We can see that the {d,p) 
and (d, n) breakup probabilities are similar in shape and 
absolute magnitude. The elastic breakup probability is 
much smaller and shows a different behavior as a function 
of £. The cross sections for the inelastic (d,p), inelastic 
(d,n) and elastic (d,pn) breakup processes are 290 mb, 
284 mb, and 122 mb respectively, which lead to a total 
deuteron breakup cross section of 698 mb for this target. 
The total (d,p) breakup cross section (which is the sum 
of inelastic (d,p) and elastic (d,pn) cross sections) is 412 
mb which is in reasonable agreement with the measured 
value of 481 mb ]3^ . The remaining difference between 
the experimental and theoretical cross sections, to some 
extent, may be attributed to the likely contributions of 
quasi-elastic channels like (d,dp) to the inclusive proton 
spectra. 

An important feature of Fig. 1 is the fact that for £ > 
21, the elastic breakup probability (EBP) exceeds the to- 
tal transmission coefficient. This shows the inadequacy 
of the phenomenological deuteron optical potential (used 
to determine the transmission coefficient) for larger val- 
ues of £. The Woods-Saxon shapes used to obtain such 
optical potentials are too restricted in their radius de- 
pendence. It is obvious that these potentials do not ac- 
count correctly for the elastic breakup channel (which we 
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shall refer to as dissociation in the following) at large dis- 
tances. We also note from this figure that both elastic 
and inelastic breakup probabilities are relatively large at 
very small partial waves. This is the consequence of weak 
absorption of neutrons and/or protons. For breakup re- 
actions involving heavier projectiles, these probabilities 
are quite small at these partial waves |12^,S3| . 
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FIG. 1. Calculated breakup probabilities in the scattering 
of 56 MeV deuteron from ''^V. The open triangles show the 
transmission coefficients calculated with the optical potentials 
as explained in the text. 

Now we proceed to determine the potential due to the 
dissociation process. As already pointed out in Ref. [p^ , 
it is the Coulomb force which is responsible for the long 
range part of the EBP. In case of pure Coulomb interac- 
tion, the eikonal approximation may be used to relate the 
imaginary phase shifts to an imaginary potential p^ , ^ . 
Since (Jb~up,d defines uniquely an imaginary phase-shift, 
we can use this method to determine the imaginary pure 
Coulomb dissociation potential by fitting the EBP cal- 
culated by the post form DWBA theory with Coulomb 
interactions only between deuteron and proton and the 
target nucleus. 

On the left side of Fig. 2, the Coulomb EBPs arc shown 
with Coulomb interactions obtained from: (1) an uni- 
form charge distribution (with a radius (i?c) of 1.3 fm) 



(solid circles) and (2) a Woods-Saxon charge distribu- 
tion with radius and diffuseness parameters of 1.3 fm and 
0.65 fm respectively (open circles). We observe that the 
Coulomb dissociation is sensitive to the charge distribu- 
tion for small i values. For larger partial waves, however, 
they produce identical results which is to be expected. It 
is explicitly clear (together with Fig. 1) that the Coulomb 
force determines the EBP at higher partial waves. We 
have used EBP calculated with the Coulomb potentials 
generated from the Woods-Saxon charge distribution to 
perform the fits. In fact, this diffuse Coulomb potential 
has been used in all the post form DWBA calculations 
presented in this paper. 




r (tm) 



FIG. 2. Results of calculations without nuclear inter- 
actions. The left part of the figure shows the calculated 
Coulomb elastic breakup probabilities for the same reaction as 
in Fig. 1 obtained with Coulomb potentials generated with a 
uniform charge distribution (solid circles) and a charge distri- 
butions with Woods-Saxon shape (open circles). The dashed 
and solid curves show the best fits obtained by using 2 and 6 
parameter parameterizations. The right part shows the corre- 
sponding best fit potentials obtained by fitting to the breakup 
probability calculated with the Woods-Saxon charge distribu- 
tion. The solid (dashed) curve represents the results obtained 
with six (two) parameter parameterizations. 

We used a two parameter, (ae~'"'), as well as a six pa- 
rameter parameterization, [(ai -f a2r"'^/^ + a3r+a4r^)/ (1 + 
air)\e^~^^^~^^^ \ in our fitting process. The best fit 
breakup probabilities obtained with the two procedures 
are shown by dotted and full curves in the left part of Fig. 
2. The corresponding imaginary potentials are shown on 
the right side of this figure. The six parameter search pro- 
cedure (with hi and &2 being equal to 0.1118 and 0.0113) 
provides a better fit to the EBP. The two parameter fits 
had a radial dependence of e^'^'^''. In comparison to 
this the Woods-Saxon potential behave approximately as 
g-r/a (^.^yj^]-^ q being the diffuseness parameter), 

at larger distances. Therefore, it is not surprising that 
the transmission coefficients calculated with conventional 
optical potentials drop too fast at larger partial waves. 

It is clear that potential due to the Coulomb disso- 
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ciation (PCD) is very weak. We found that the elas- 
tic scattering angular distributions calculated with phe- 
nomenological optical potential were almost unaffected 
by inclusion of the PCD to its imaginary part. Also this 
leads to a change of less than 1% in the total reaction 
cross section. The transmission coefficients at lower par- 
tial waves are affected by the inclusion of PCD, but the 
nuclear effects damp them out strongly in this region, as 
is discussed in the context of Fig. 3. 

We, therefore, added the PCD to the phenomenologi- 
cal imaginary potential as it leads to the total transmis- 
sion coefficient which is larger than (or equal to the the 
EBP) for all the partial waves as can be seen from the 
left part of Fig. 3, where the total transmission coeffi- 
cients obtained with (solid triangles) and without (open 
triangles) adding the PCD to the phenomenological op- 
tical potential are shown, ft can be seen that, addition 
of PCD to the phenomenological potential removes the 
anomaly described above. In the right side of Fig. 3, we 
compare the results of the Coulomb EBP obtained with 
(by using Eq. (13)) and without nuclear distortion. We 
note that nuclear distortion effects suppress the Coulomb 
EBP strongly for lower partial waves. However, for larger 
partial waves the two calculations produce identical re- 
sults. 



i 

+exp{-{{r + R,)hY)] (3.1) 

The value of 7 was taken to be ^2/3 [|6| , and a total of 
15 terms were included in the sum. The coefficients Ai 
and positions Ri were varied in order to get the best fits 
to the elastic, inelastic and total breakup probabilities, 
which are shown, in the left, middle and right parts of 
Fig. 4 respectively. The solid circles show the results of 
post form DWBA (same as that in Fig. 1), while open 
circles represent our fits. In case of EBP, we also show 
the results for pure Coulomb case (open triangles) (the 
same as shown in the right side of Fig. 2). We can see 
that Coulomb EBP agree very well with that obtained 
by fitting the total EBP for values of £ > 21. It should 
be remarked here that in the study of Christley et al , 
who have used a similar method to determine the fusion 
potential by fitting to the fusion cross sections calculated 
within a coupled-reaction-channel model, the quality of 
fits were not as good. They attribute this failure to the 
lack of ^-dependence in their fitting potential (which were 
taken to be of the Woods-Saxon type). However, we do 
not require such a dependence in our fitting procedure. 
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FIG. 3. The left part shows the transmission coefficients 
in the elastic scattering calculated with the phenomenologi- 
cal optical potential with (solid triangles) and without (open 
triangles) addition of the imaginary potential due to the elas- 
tic Coulomb dissociation. The right part shows the breakup 
probabilities for elastic Coulomb dissociation calculated in the 
eikonal approximation (same curve as shown in figure 2) (solid 
circles) and those obtained in the full calculations with nu- 
clear distortions included (open circles). 

In order to find the potentials due to the inelastic 
breakup and nuclear dissociation, we used Eq. (13). 
The inelastic and elastic breakup probabilities calculated 
within post form DWBA (Fig. 1) were fitted by using a 
sum of the Gaussians (SOG) parameterization Q 
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FIG. 4. The elastic, inelastic and total breakup prob- 
abilities calculated within the post form distorted-wave 
Born-approximation (solid circles) and our best fit values to 
them (open circles). Open triangles, in the left part of the fig- 
ure, represent the breakup probability for the elastic Coulomb 
dissociation (same as those shown in the right part of fig. 3 
by open circles). 

The potentials giving these breakup-probabilities are 
shown in the upper part of Fig. 5. The solid, dashed 
and dotted curves show the potentials due to the inelas- 
tic breakup (which has been plotted after multiplying the 
actual values by 0.09), nuclear dissociation and Coulomb 
dissociation (already shown in Fig. 2) respectively. We 
note that the inelastic breakup and nuclear dissociation 
potentials are strongly concentrated in a region around 
6.5 fm and 8 fm, respectively. Both potentials have about 
the same width, 2 fm (FWHM), which is very similar to 
the rms-radius of the deuteron. The difference in the lo- 



6 



calization of the inelastic and dissociation potentials con- 
firms the fact that inelastic breakup of the projectile oc- 
curs in regions closer to the target nucleus as compared to 
the elastic breakup. It is worthwhile to note that the po- 
tential due to the inelastic breakup is quite large. There- 
fore, those calculations [|o|jl^j3^,Q where this mode of 
breakup is not considered are likely to produce weaker 
imaginary potentials due to the breakup channels as com- 
pared to ours. Due to its long range, the potential due 
to the Coulomb dissociation is still of appreciable magni- 
tude at large distances. Coulomb breakup thus, can take 
place even outside the charge distribution of the target 
nucleus. 

In the middle part of Fig. 5, we compare the imaginary 
part of the phenomenological optical potential (solid line) 
with the bare potential (obtained by subtracting the po- 
tential due to breakup from it) (dashed line). It is clear 
that breakup is the dominant absorption effect in the 
surface region. The bare potential could be the starting 
point of double folding model calculations. 
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r (fm) 

FIG. 5. The upper figure shows the imaginary potentials 
due to inelastic breakup (arbitrarily multiplied by 0.09 to fit 
in the plot) (solid line) and the potential due to nuclear dis- 
sociation (dashed line). The dotted curve shows the Coulomb 
dissociation potential. The middle part shows the optical po- 
tential with (solid curve) and without (dotted curve) effects 
of breakup. In the lower part the dashed, dotted and solid 
curves show the neutron, proton and matter densities for ^^V 
respectively, as calculated by Fayans E^. 



In the lower part of Fig. 5, we show the radial depen- 
dence of the neutron, proton and matter densities for the 
^^V nucleus WE note that breakup probabilities 

peak in a region where the matter and charge densities 
are much smaller as compared to the central density of 
the target nucleus. Therefore, breakup process is really 
sensitive to the extreme peripheral regions. This also 
confirms the spectator role of the target nucleus in the 
breakup process (Serber picture). However, the region 
where the breakup potentials peak is much larger than 
the half density radius deduced by Serber |Q . 

It should be stressed that our method of obtaining the 
imaginary part of the optical potential due to the breakup 
channels is phenomenological in nature. However, it has 
several good features. The distinctly different nature of 
the total breakup cross section as a function of the in- 
cident partial waves is automatically taken into account. 
The Coulomb breakup process which is responsible for 
the long range part of the elastic breakup probability 
(and the dissociation potential) is included into our cal- 
culations. Furthermore, we include the inelastic breakup 
mode which makes the largest contribution to the imagi- 
nary potential. Both these effects were ignored in calcu- 
lations presented in Refs. pO| , pT| . 

IV. SUMMARY AND CONCLUSIONS 

In this paper we calculated the imaginary part of the 
optical potential due to the breakup channels by using 
a phenomenological method in which contributions to 
the reaction cross sections from the elastic and inelas- 
tic breakup cross sections of the projectile in the field 
of target nuclei, calculated in the post form distorted- 
wave Born-approximation, are used as constraints in a 
fitting procedure. In this method, only the first order 
breakup process (the coupling of the elastic channel to 
the breakup channel) is considered. Then unitarity of the 
^-matrix is used to determine the influence of breakup 
on the elastic channel. 

The Coulomb part of the imaginary dissociation po- 
tential (obtained from the pure Coulomb elastic breakup 
probabilities) accounts for the long range part of the ab- 
sorption and removes the anomaly where for the partial 
waves beyond 20, the elastic breakup probabilities were 
found to be even larger than the transmission coefficients 
calculated with the usual phenomenological optical po- 
tentials. To our knowledge this is the first calculation 
where the absorption due to the Coulomb dissociation 
process has been included in the optical model. The po- 
tentials due to the nuclear dissociation are found to be 
peaked in the region around 8.0 fm. 

The imaginary potential due to inelastic breakup is 
about an order of magnitude larger than that due to 
the dissociation process, and is concentrated at some- 
what shorter distance (around 6.5 fm) as compared to 
the latter. This suggests that the inelastic breakup takes 
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place at distances closer to the target nucleus in compar- 
ison to the elastic breakup. The magnitude of the total 
imaginary potential, due to the breakup channels (IPBC) 
is, therefore, almost solely due to the inelastic breakup 
process in the region around 6-8 fm. This is consistent 
with the dominant contribution of this mode to the to- 
tal breakup cross section. This is the main reason for 
our IPBC being stronger than the imaginary part of the 
dynamical polarization potential due to breakup chan- 
nels calculated by other authors, who have ignored the 
inelastic breakup mode. 

Our method can be applied to any projectile having 
a strong breakup channel. Therefore, It would be in- 
teresting to use this method to the elastic scattering of 
halo nuclei (eg. ^^Li and ^^Be) where the breakup 

cross sections are significantly enhanced and the effect of 
breakup on the elastic scattering is very strong. 
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